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Abstract 



A proposition wherein the dynamics of dust during mini-rover operations on lunar surface taking into account 
of probable resistance offered by an equilibrium local concentration of gases of near-surface lunar atmosphere 
can be best described by considering a weak linear drag force in the equations of motion is presented. The 
equations of motion vis-a-vis the one considering a vacuum near-surface lunar atmosphere are transcendental 
in nature whose exact closed-form solutions can be obtained using Lambert W-function only. An analysis is 
made giving specific consideration of Chandrayaan-2 Rover and results are presented. 
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1. INTRODUCTION 

The lunar surface is known to be covered by a thick 
blanket of fine dust. The lunar “regolith” is a blanket 
of soil and rocks 3 to 20 m thick that covers the lunar 
surface. The size of lunar soil particles varies from 
50 pm to 1 mm in diameter[l]. Through various lunar 
soil mechanics experiments, it is reported that lunar 
soil exhibits strong cohesive forces. The predominant 
forces between lunar particles are van der Waals and 
electrostatic interactions. The lunar dust consists of 
extremely fine and electrostatically charged particles. 
Lunar particle cohesion has been estimated to be 
roughly 0.035 to 0.05 N/cm^ from lunar surface 
experiments [2]. Lunar dust clings tenaciously to 
surfaces and during Apollo missions, adhesive forces, 
were observed to vary roughly from 0.01 to 0.1 
N/cm^ and deposition of such dust on sensitive surfaces 
like solar panel, optical lenses etc. causes serious 
problems [3]. In this regard, an analysis is carried out 
to understand the dynamics of dust under the 
equilibrium near-surface lunar atmosphere during mini 
rover operations giving specific consideration of 
Chandrayaan-2 rover and results are reported. 

2. ROVER OPERATIONS AND ROVER 
WHEEL MECHANICS 

During the rover operations on the surface of Moon, 
dust will get collected on the rover wheel and 
subsequently, it will either remain stuck to the wheel 


or be kicked off following the trajectory given by the 
dynamics of rover wheel depending on the mode of 
rover operation. The modelled analysis of the dynamics 
of dust can be studied under the following cases: 

i. Normal rover operation (i.e. rolling wheel); & 

ii. Contingency rover operation (i.e. slipping 
wheel). 

2.1 NORMAL ROVER OPERATION 

During the normal mode of rover operation, rover 
wheels roll on the surface. Motion of a rolling wheel 
can be considered as the combination of pure rotational 
motion of the wheel and pure translation of the wheel 
about the centre of the rover wheel. Let us consider 
a rover wheel of radius “R” rotating with an angular 
velocity of “co” as shown in Fig. 1. Let the velocity 
of the centre of the wheel be Then, the 

condition of pure rolling gives 

Vcm=coR ( 1 ) 



Fig. 1 Fig. 2 
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The velocity of any point “P” on the circumference of 
the wheel is different from and can be calculated 
as follows as shown in fig. 2. A rolling point “P” is 
acted by two component velocities; one component 
velocity due to rotation and other component 
velocity V due to translation. Thus, the velocity of 
the point “P” can be written asFp = Fc + Pure rolling 
condition yields: ||71 = = \Ycm \ ^\^c\ 

= ( 2 ) 

Thus, a dust particle collected at “P” attains the 
velocity given by (2) and will tend to fly tangentially 
along a trajectory which will be subsequently 
shown to follow the equation (10). To study the 
trajectory of the dust, the resultant velocity given 
by equation (2) can be resolved into two orthogonal 
axes defined as shown in fig. 1. Thus, it can be written 
as 


Vp = + Vpp ; where 

Vpx =u^= VcM (■v)+ VcM cos6> (- x)= VcM (l - cose)i 


VpY ~^y ~ ^CM sill® (y) 


( 3 ) 


Now, to find the direction of the resultant velocity at 
“P”, let “a” be the angle made by the resultant velocity 
“Vp” with the horizontal direction. Then, 


tana 


Pr 1 in 6 

= Y-l = cot - = tan- 


u 


12 


= -(7r-e) 
2^ ^ 


i.e. resultant velocity at “P” bisects always the angle 
in -0) . 


2.2 CONTINGENCY ROVER OPERATION 

During this rover operation mode, wheel(s) of rover 
is (are) assumed to get stuck on (dusty) ground so 
that rover wheel is in slip operation mode. In other 
words, rover wheel is experiencing only 100% rotation 
and 0% translation. Since there is no translational 
motion of the wheel, the velocity of a point “P” is 
given by 

^Vp =V^ =GyR (4a) 

Thus, unlike a rolling wheel, all the circumferential 
points on a slipping wheel will have the same velocity 
whose magnitude is given by equation (4a). Thus, a 
dust particle adhered at P will tend to fly tangentially 
with initial velocity given by (4a) along a trajectory 
as derived in equation (10). 


The components of the resultant velocity obtained in 
equation (4a), resolving along the two orthogonal axes 
as defined in fig. 1, are given by 

VpY = (oRsinQ (+ j)) and Vpj^ =Ux = coi^cos0(- jc) 

(4b) 


3. DUST DYNAMICS 

It is a well known fact that the acceleration due to 
gravity on the surface of the Moon is about one-sixth 
of that on Earth: roughly 1.63 m/s^. This is too weak 
to retain an atmosphere like that of Earth. However, 
it is capable of retaining a thin atmosphere which 
consists of traces of Ne, He, H^, and Ar, and even 
smaller amounts of CO^ and H^O. Some of these gases 
are trapped in the soil during the cold lunar night and 
escape as the surface heats up from the sun. Gas 
abundances of all the constituents of lunar atmosphere 
and some scale heights in the lunar atmosphere as 
well as the variations in abundances of the four 
principal lunar atmospheric gases (Ne, He, H^, and 
Ar) over day/night cycles at different latitudes and 
longitudes are available [4]. The daytime undisturbed 
gas concentrations (molecules/cm^) of the four major 
constituents of lunar atmosphere viz. ^^Ne, He, H^, 
and "^^Ar are respectively -4x10^ to 10"^ ; -8x10^ to 
4.7x10^; -2.5-9.9x10^; and -2x10^ [4]. The respective 
numbers during night time are 10^; 4 to 7x10"^; 10"^ 
to 1.5x10^; and <10^ [4]. Also, density of gases 
decreases with increase in altitude. Thus, these local 
concentrations of trace gases may be assumed to offer 
resistance to the motion of minute lunar dust particles 
during mini-rover operations like Chandrayaan-2 rover 
of height less than a meter. Thus, the dynamics of 
dust during rover operations on lunar surface taking 
into account of probable resistance offered by an 
equilibrium local concentration of gases of near-surface 
lunar atmosphere might be best described by 
considering a weak drag force in the equations of 
motion.. The weak drag force can be assumed as a 
drag force which is linear speed and its general form 
can be written as = -bv, where “b” is a constant 
whose value depends on the properties of the 
medium and on the shape & dimensions of the 
object. Eor a spherical object, the constant “b” takes 
Stoke’s constant 67ir|jr, where r is the radius of the 
sphere and is the dynamic viscosity of the fluid 
(medium). 
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To determine whether linear or quadratic drag applies 
to a system, one must calculate the Reynold’s 
number, Such a linear (in speed) drag force actually 
occurs for low Reynold’s number, 51g<l and 
represents a streamline flow. For l«9f^ «10^ , the 
motion is turbulent and quadratic drag applies. 
Quadratic drag is suitable for motion of dust particles 
in Earth’s atmosphere with thick atmosphere and also 
in Mars atmosphere wherein the dynamics of dust is 
very much affected by the wind velocity which varies 
from 2 m/s to 10 m/s near the surface [5]. The motion 
of lunar dust, on the other hand, is mainly gravity 
driven. Though lunar dust are available in different 
textures, lunar dust particles are assumed to be 
spherical for simplicity and average lunar particle 
density is taken to be 3.01 gm/cm^[3]. Then, for a 
spherical dust of radius r moving with speed v in the 
local cloud of gases (fluid) of density and dynamic 
viscosity (or kinematic viscosity r]^), the Reynold’s 
number and drag force are given by 




_ 2rv _ 2rpfV 
- ^Tirj^rv for 


m,<i 


(5a) 

(5b) 


The terminal speed is then given by 

4 3 

6nrijrVf = mg„, = -nr p,g„, (5c) 

where p^ is the density of a spherical dust and is 
the lunar gravity. In order to apply equation (5 a) to 
the whole of the trajectory, we require 91 ^ < If or speeds 
up to the terminal speed. For solving equation (5c) for 
speed up to terminal speed i.e. v^v, the radius of 
the spherical lunar dust is found to be 


4p fPsgM 

Assuming motion of a lunar dust in a local 
concentration of gases of each of the major constituents 
of near-surface lunar atmosphere, for^I^ < 1, one finds 
that equation (5d) applies when 


Radius (pm) Ne 

He 

Ar 

r < 8.18 

< 101.8 

< 51.6 


Spheres this small is of our interest as scale of lunar 
dust size is of the order of these spheres. Thus, linearity 
of drag force may be considered in the equations of 
motion. 


(5d) 


3.1 VECTOR EQUATIONS OF MOTION 

Defining the +Y-axis along the upward vertical 
direction, and +X-axis along the horizontal right 
direction as shown in fig.l, the equations of motion 
of a dust particle of mass “m” at point “P” can be 
written as follows: 

dv^ 7 

= and m — - = -mg-bv (6) 

dt dt 

The first order (in velocity) uncoupled linear 
differential equations (6) can be effectively combined 
in a single vector equation (7a) as follows: 

m^=F + l (7a) 

dt 

Where the gravity force F and drag force f^ are given 
by 

F - -mgy and fd - -bv (7b) 

(7c) 

dt 

This non-homogeneous first order differential equation 
is the desired equation whose solution is sought. And 
it can be solved as follows: 




(7d) 


Change of variable: u=v -reduces equation (7d) 

b 

to: 


du , ^ 
m — = -bu 
dt 

The solution of equation (7e) is given by 


^ F - —t 
^ V = —Ae ^ 
b 


(7e) 

u-Ae^ 

(7f) 


Evaluation^ using the initial condition i.e. t = 0; v = u, 
equation (If) becomes 


^ cff(t) F 

^V = —^ = —+ 

dt b 


F m 

^{tU—t - 

b b 


^ F 

u - 

b 




/ 




+ B 


( 8 ) 


(9a) 


Evaluating again B using initial condition i.e. t = 0; 
r = 0, the vector solution of equation (7d) in r is 
obtained 


b b 


f 




U-- 


\-e 


t 


(10a) 


3 























Since F = scalar equations of equation (10a) are 
given by 


m 


b 


\-e 


b \ 
- 1 


b b\^b 


\-e 


(10b) 


(10c) 


If “u” is the initial velocity with which a dust particle 
leaves the rover wheel making an angle “0” with the 
horizon, then the dust particle will traverse a trajectory, 
whose parametric equations are given by equations 
(10b) and (10c). From equation (10b), the time “t” is 
obtained as 


m , 

t = -In 

b 


f 


1 - 


bx 


V 


(11a) 


Substituting equation (11a) in equation (10c), the 
equation of complete trajectory of dust particles is 
given by 




bx 


(lib) 


Since we get now the parametric equations (10b) and 
(10c), the dynamics of dust particles can be studied in 
terms of its time of ascend, maximum height, time of 
flight, range etc. 

3.2 TIME OF ASCEND 

Mathematically, time of ascend is defined as the time 
for whichy(/)= 0. From equation (10c), we get 

The condition 0 yields the time of ascend t^ as 

follows: 

b 

- h 




m. \ ^ b 

h=—^ 1 +-Wy 

b \ mg 


3.3 MAXIMUM HEIGHT 


( 12 ) 


Mathematically, maximum height H is obtained from 
the condition: H = Thus, from equation (10c), 
we get 

mg m( bng \ ^ 

^ H — - U H-H-1 — C 

b ^ b[^ b 


b \ 


(13a) 


Substituting the value of t^ from equation (12) in 
equation (13a), we obtain the maximum height H 
attained by dust particles as follows: 


f m 

f. b \ 



U, 

j ^ mg ^) 


m mg 
b[ ^ b 


1 - exp<^ - In 


1 ^ 

1 - bly 

mg 


(13b) 


3.4 TIME OF FLIGHT 


Mathematically, time of flight is defined as the time 
for which y{t) - 0. Thus, from equation (10c), we get 


mg m , 


I bng 

b 

g 


mg 
^ b 


l-e 


^2 +1 


= 0 


(14) 


Equation (14) is an implicit equation whose solution 
cannot be expressed in terms of elementary functions. 
One has to solve this transcendental equation using 
numerical techniques or approximations. However, our 
interest is to solve equation (14) analytically. This is 
possible only with Lambert IT-function. This paper 
presents an approach to use Lambert IT-function 
for rover applications and it describes finding the 
exact closed-form explicit solutions of probable 
transcendental equations which govern the dynamics 
of dust in the near-surface trace gases lunar atmosphere 
during lunar mini-rover operations, using IT-function. 

3.5 LAMBERT W-FUNCTION AND ITS 
OVERVIEW 

Lambert W-function commonly known as “W-function” 
originated from the work of J. H. Lambert [6] on 
trinomial equation that he has published in 1758 and 
was discussed by Euler [7] in 1779. Although rarely 
used, its properties and algorithms are well documented 
[8-12]. Nowadays, it is included as an in-built library 
function in many computing systems such as 
LambertW[x] in Maple® and ProductLog [x] in 
Mathematica®. Lambert W-function, represented by 
W(z), is defined as the inverse of the function 


f(z) = ze^ satisfying W(z)e^^^^ = z 
That is: Y = Xe^ ^ X = W(Y) 


(15a) 

(15b) 




















The Lambert TT-function as defined by equation (15) 
is generally complex-valued. In general, for a given 
complex number z, equation (15) will have infinite 
number of solutions and is therefore multivalued 
function. Thus, in general solutions of equation (15) 
are written as (z) where the integers k = 0, 
±1, ±2, ... denote different branches of the Lambert 
W-function. 

For our purpose, we are interested in the real-valued 
W(x). Thus, when z is a real number x, equation (15) 
can have either two real solutions Wq(x) and W ^(x) 
for -1 / ^ < X < 0, one real solution Wjix) for x > 0, 
or no real solutions for x<-l/e. By convention, the 
branch satisfying W(x) > -1 is denoted by Wq(x) and 
is defined as the principal branch while the secondary 
real branch satisfying W(x) < -1 is denoted by 
IT ^(x). The two real branches are shown in fig. 3. 



Fig. 3: The two real branches of Wjj(x). The solid line 
shows Wjj(x) while the dashed line shows W ^(x). [11] 


Some of the values of W-function are given in the 
following table 1. 


Table 1. Some of the values for the Lambert 
W-function. Of the infinite number of branches, 
only two branches are tabulated. 


X 

W„(x) 

Wj(x) 

e 

1 

-0.5321 - 4.597i 

1 

0.5671 

-1.534 - 4.375i 

-1/e 

-1 

-1 

-1/4 

-0.3574 

-2.153 


3.6 GALILEAN TRANSFORMATION 

In order to make use of Lambert W-function for our 
specific purpose, let’s look for the general solution of 
an implicit equation given by 

pax+b = cx + d (16) 

If equation (16) can be expressed in the form tp^ = R, 
then solution involving W-function is obtained as 
follows: 


tp^ = R (17a) 

^ Inp = Rlnp 


^ t = - Lambert 

In p 


w{R\np) 


(17b) 


Without the loss of generality, to express equation 
(16) to the form of equation (17a), the following 
Galilean transformation is made. 


t = mx + n 

Thus, equation (16) becomes: 


(18a) 


. — =tp~"' +(d-—]p~'” — 

c [ m J c 


(18b) 


Equation (18b) will reduce to equation (17a) if 

A , nc ^ md 

m = -a and d -= 0; ^n = - 

m c 

Thus, the required Galilean transformation is 


t = -ax - - 


ad 


(18c) 


Hence, equation (18b) is reduced to the following form 


tp^ = R where R = -—p ^ 

c 


The solution of this equation is given by equation 
(17b) i.e. 


-ax- 


ad 


1 

In p 


Lambert W 


, ad 

a ^ - 

— p ^ \np 
c 


X = 


a\np 


- Lambert W 


a\np 
- -P 


ad 


d 

c 


(18d) 


Equation (18d) is the exact and explicit solution of 
the implicit equation (16). Thus, the exact and 
explicit form of time of flight of dust particles 
as given by equation (14) in terms of Lambert 
W-function can be written, using equation (18d), as 
follows: 
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= — Lambert W 
^ b 




1 ^ 

1 h - u 


1 ( mg 


Lambert W (x)+ — 


Where x = -| l + -^i/ expj- l + -^i/ 

m J \ V m 


3.7 HORIZONTAL RANGE 


(19) 


4.1 NORMAL ROVER OPERATION 

As derived in equation (3), in this mode of rover 
operations, we have the following component 
velocities: 


u = (1 - cos 6 ) and u = sin 6 

X CM ^ '' V CM 


(22a) 


Thus, substituting equation (22a) in equations (12), 
(13b), (19), and (21), we obtain the following set of 
equations which govern the dynamics of dust particles 
kicked off the rover wheel: 


Horizontal range R can be evaluated using either of 
the following definitions: 

R = x{t^\ or range R is the value of x for which 
y{x) = 0. 

It is convenient to determine range R using the second 
definition i.e. using the trajectory equation (11b). 
Therefore, equation (11b) reduces to 


= -^ln[ 1 + — VcM sin0 
b i mg 


= Jl + A(/c^sine 


m { 

~b 


(22b) 


flcMsinO +-^1 

1 - exp<^ 

-ln[l + — sine ]i 

1 , b J 

1 

[[ mg )\\ 


(22c) 


If mg] m^g. L bx 

— bc + —^In 1-- 


mu^ 


= 0 


. , 1 
‘ -- 

m gu^ 




mu^ 


( 20 ) 


Equation (20) is again a transcendental equation whose 
explicit solution can be written in terms of Lambert 
W-function only. Thus, the horizontal range R is 
obtained as follows using equation (18d) 


R = - 


, Lambert W 





f buA 

( bu.. ^ 

■ 

- 1 + ^ 


_1 



mu^ 


R = Lambert W(x)+ 

^ mg ] b 


( 21 ) 




Where ^ - 


1 ^ 

1 H- U 

, mg J 


' J ' 

exp^- 


mg )\ J 


4. CONTEXT OF ROVER OPERATIONS 

Having arrived at the general formalism of dust 
dynamics, now let’s apply it to study the dynamics of 
dust specific to the following modes of rover 
operations as discussed in preceding sections. 


t^ =—Lambert W 
^ b 


- 1 + 


bVcM sin0 
mg 


exp<- 1 + 


sine 
mg 


+ Lv^^sinO + ^ 
gK b 


(22d) 




Z)2fFcMSin0+^l 


- 1 + - 


^'exp 


bVcM sin0 


mg 


I I ^IcM sin0 
mg 


+ fl"cM(l-cose) (22e) 

4.2 CONTINGENCY ROVER OPERATION 

In this mode of rover operation, as derived in 
equation (5), we have the following component 
velocities: 


u = (dR COS0 and u = coR sin0 


(23a) 


Thus, substituting equation (23a) in equations (12), 
(13b), (19), and (21), we obtain the following set of 
equations which govern the dynamics of dust particles 
kicked off the rover wheel: 


= —In 1 + -^coRsinO | 
b [ mg J 


(23b) 
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(—1 lnfl + -^o_Rsin0l 

^ h ~ hvac 

1 “ :rf~\ivac + ^'(sh 

\b) y mg ) 


2ym) 


(24b) 


m ( 

J 


[co7?sin0 

1 - exp^ 

- In [ 1 + -^ CO 7? sin 0 11 

1 b J 

1 

1 1 'wg Jj_ 


(23c) 


Thus, we see from equation (24b) that trace gases in 
lunar atmosphere decrease the time of ascend of dust 
particles as expected. 


^9 = —LamhertW 


I, bcoRsinO] I T bayRsinQ ] 
-I 1 +- exp^- 1 + 


mg ) 


mg ) 


+ —[ coi^sind + 


R = _P^^^S^^Lambert W 
b^LiRsine +'yj 


(23d) 


Substituting £ -- u in the equation (13b) for maximum 

mg 

height H and expanding in Taylor series form, we get 


r ^ ^coi^sinG ^ 


1 - 1 1 - In (l + 8 )+ — ln(l + 8 )ln(l + 8 )-... 

^1] 

1 mg j 





exp<- i + 


Z)co7^sin0 j 
mg J 




£ -+ .. 

2 , 


m 

+ —coi^cosG 
b 


(23e) 


5. EFFECTS OF LUNAR TRACE-GASES 
ATMOSPHERE 

As discussed above that the non-zero gravity of the 
Moon is manifested to a non-vacuum atmosphere 
consisting of trace gases. In order to see the effects of 
these local trace gases near the lunar surface, we have 
to see the Taylor series form of equations (12), (13b), 
(19), and (21). We have the time of ascend t^ as given 
by equation (12): 


m, 

U = —In 
^ b 


1 ^ 

1h- u 

mg J 


^ m. ^ \ . b 


= —ln(l + e) where e =- u and 

b mg ^ 


b < 10“^ gm / sec for a lunar particle of diameter < 1 mm. 
Using Taylor series expansion, we get 


3 4 

8 8 8 

8 -+-+ . 

2 3 4 


m 

^ u = —s 


- — 


1 - 


8 8^ 8^ 




(24a) 


H = 


H=£ 


mUy \ b m 1 bUy m mg 

b 2 mg ^ b ^ 2 mg b b 

\ 2 m( mg^^ 

+ -£ — +.... 

4 b[ y b 


1 mu \u 1 2 ^ f mg 

-^-^ + -£ —\u^+—\ + . 

2 b 2 g 4 b{^ b 


1 


1 mUy \ Uy 1 2 ^ f mg 


H- £ — \ U., -+_ 


lb IgA b 

^_Uy b Uy Uy ^ I Uy ^ I b ul l ^ 

2g m g 2g 4 4 m g^ 4 m g^ 


1 K r.^r.4\ 

+ 7 —^ + 

4 m g 


H = ^ 
2g 


H = W, 


H = H,. 


f„ \ 


, bu 

1 -+ — 

mg m 


1_A 2 

^ hvac 2 ^ 

m m 



b 

^ hvac 
m 


(25a) 


(25b) 


Thus, it is seen from equation (25) that presence of 
trace gases in lunar atmosphere decreases the maximum 
height attained by dust particles. 


Since / = ^ is the time of ascend of dust particles in 

4vac ^ 




vacuum, equation (24a) reduces to 


Total time of flight t^ is obtained from equation (19). 
One can expand the W-function in series using the 
series expansion of W-function. However, just to see 
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the first order effect, the Taylor series form of time of 
flight as given by equation (14) may be used. 

^ 1 - — ^2 +--^^2 +^(^ 2 )= “7 - - - +1 


m 


m 


mg 


g 


1 

2m^ I mg 

^ b 


m 


= 0 


(26a) 


Hence, we see from equation (27b) that horizontal 
range of the dust particles decreases as expected. 

Thus, presence of trace gases in the near-surface lunar 
atmosphere due to non-zero gravity of the Moon 
modifies the dynamics of dust during rover operations 
and the worst case of the dynamics is given by the 
ideal vacuum conditions of Moon i.e. considering as 
if Moon does not have an atmosphere at all. Thus, 
following results are applicable. 


The non-trivial solution of equation (26a) is given by 


2m^ 

b 

g 

2i/ 

y 

1 

b^ 

m 


g 

1-1 

+ 


‘-Ivac 


1 b 

1 h —61 


m 


(26b) 


Thus, we see from equation (26b) that the total time 
of flight of dust particles also decreases because of 
the presence of drag force due to lunar trace gases. 

Again, using the same analogy as in the time of 
flight, we take Taylor series expansion of the 
equation (20). 


^1- 


1 

m^g 




1 b‘' 


b ) 2 m^g^u 


4„2_2 y 


mg^ 

' b ] 




-R + \ 


mUr 


b^ R ( bR 1 b^ 


m g 

b^ 


u,, +- 


o 4 2 2 I 




m^gu^ 


R 


^ 2m^gu 


I - ^Uy+mgjR 


= 0 


(27a) 


The non-trivial solution of equation (27a) is given by 

^ -Uy + - J - ^Uy + mgj R = 0 

2m gu^ 


R = 


2m^gu^Uy 

2m^gu^Uy 

2U^Uy 

1 

]^Uy + mg] 

2 2 
m g 

+ 

^ g 

1- 

+ 


mg _ 


m g _ 




R. 


\ b 

m 


n2 


(27b) 


The maximum height attained by dust particles during rover 

2 


operations: ^max - 


u 


y 


2g 


The maximum range of dust particles: 

2Uy.u,, 

D _ 

^max 

g 

The total time of flight of dust particles: 
2u 

T = —- 

^ max 


(28a) 


(28b) 


(28c) 


g 


6. MINIMUM VELOCITY OF WHEEL TO 
BREAK A DUST PARTICLE OFF THE 
WHEEL 

Assuming that a dust particle of mass m and radius 
r sticks to the rover wheel of mass M and radius R at 
the point of contact A and it is transported at the point 
P traversing a distance arc AP during time t. Also, 
let’s assume that dust particles adhere to the surface 
of wheel with a force ‘F’ and the co-efficient of friction 
between the wheel and dust particles is p. Let I be 
moment of inertia of wheel, F be the radial force and 
N be the normal reaction. 


If V and CD be the linear velocity and angular velocity 
of centre of mass (CM) of the wheel before picking 
up a dust particle. After a dust particle gets stuck to 
the surface of wheel, angular velocity of the centre of 
mass will be changed to CD^ and V becomes V^. In 
ideal situations, rolling motion of wheel on undulating 
dusty surface will cause continuous variation of CM 
in vertical direction. However, for simplicity in the 
calculation, this is neglected and the position of CM 
of wheel is assumed to be at a constant height from 
the ground. Frictional force of the ground surface will 
try to cease the rolling motion of the wheel after some 
time. However, the force of traction compensates this 
frictional force and maintains uniform angular velocity 
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of the CM. Moreover, since we are interested in the 
dynamics of dust, we consider only the forces acting 
on the dust at the instantaneous position P, and we 
neglect such kind of forces. 




Defining the axes as shown in figure 4b, we can write 
down the equations of motion of dust particle at P as 
follows: 


7 - A. 

^ ^ =F-N-mgcosO (29a) 

ma^ = mgsinO + / (29b) 


where is the tangential acceleration of dust particle 
and / is the frictional force between wheel and dust 
particle. 


Applying work-energy theorem, we get 


Kf = -m*2V^^(l-cosQ)+-MV^^ +-l\ 


2 R 


(33b) 


Using equation (31) and (33), equation (30) becomes. 




2/77(1 - cos 0 )+ M H —Y 




= U 


M + - 




- 2mgR(l - COS 0 )- 2 iiNR 0 


(34a) 


Substituting the value of from equation (29a), 
equation (34a) becomes 


(R + r)F n(R + r) 


- g(R + r)cos 0 


2/77(1 - COS 0 )+ M + - 


+ 2mgR(i - COS 0 )+ 2\jiNR0 


M + -!— 


(35a) 


At the point of BREAK-OFF, N = 0; thus equation 
(35a) reduces to 


v = 


{R + ry 


- g{R + r)cos0 


2m(l-cos0)+M + — +2/77g/^(l-cos0) 


M + - 


(35b) 


Work done by all forces = Change in kinetic energy 
W, + W, + W + W,, + W = K - K. (30) 

The non-zero work done is obtained as follows: 


Equation (35b) defines the required minimum 
velocity “V” of rover wheel, to break a dust particle 
off the wheel at an angle “0” during normal rover 
operations. 


Work done by gravity acting on mass m 

0 0 

Wmg = ^mg.ds = j/77gsin0di(-7?0)=-/77gi?(l-cos0) (31a) 
0 0 

Work done by friction f 

0 0 

Wf = jAijV.cfi = uNR COS7C ^e=-nNRe (31b) 

0 0 

Initial Kinetic energy K 

K,=-m\v + wy.Rf +-MV^ + -I(o'^ (32a) 

2 I I 2 2 

Final Kinetic energy K^ 

Kj +©2 X 7?| +— MVy + ~ l(x>i (32b) 

For pure rolling, V= coR &V^ = co^R = thus equation 
(32) becomes 

1 / V 1 7 1 7 1 7 1 

K. =-m(V-coRj +-MV^ +-1(0^ =-MV^ +-1— 

^2^ ^2 2 2 2 (33a) 


7. RESULTS AND DISCUSSIONS 

Having derived the complete set of equations from 
the proposition which governs the dynamics of dust 
during rover operations, it can be now applied to four- 
wheeled Chandrayaan-2 rover of 18 kg mass (on the 
Earth). 

The rover is planned to land on a vast plain area near 
lunar South pole at 88°S latitude. The lunar surface 
temperature in and around the landing site latitude 
varies between -120°C to -160°C [13-17]. The rover 
has a solar panel at a height of minimum of 75 cm 
from ground and rover is configured such that 
VcM = for all modes of operation. Considering 
diameter of a dust particle as 50 pm and corresponding 
mass of a dust particle as 1 pg [3], the calculations 
are carried out. Since the properties of lunar near¬ 
surface atmosphere are not known, calculation is 
carried out for each gas component of atmosphere 
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whose property is nearly found. Though the dynamic 
viscosity of gases varies with temperature and 
pressure, its variation with pressure is very weak and 
hence is not considered. However, the variation of 
with temperature is calculated using Sutherland’s 
model by which varies with absolute temperature T as 
follows: 

3 

T„ + s f r V 

T+ s\t^) (36) 

in which Tj^is the viscosity at the reference temperature 
Tq and S is the Sutherland constant. Thus, for helium 
and hydrogen gases of known Sutherland constants, 
values are calculated using equation (36) at the landing 
site temperature of -120°C. However, argon gas of 
unknown Sutherland constant, available value at 
-100°C is being used [18]. Though Chandrayaan-2 


rover is designed for a nominal speed of 1 cm/sec, 
analysis is carried out for three different speeds and 
results are tabulated in the following tables (Table 2, 
3, and 4). 

Also, it can be calculated that speed of rover required 
to throw dust particles up the solar panel height >155 
cm/s. Again, considering a dust particle of diameter 
1mm, mass Img, and rover of each wheel mass 500 
gm, diameter 180 mm, it can be shown that minimum 
speed required to break a dust particle off the rover 
wheel > 5 m/s. 

Thus, a rover moving with a speed of 1 cm/s can not 
break dust off the wheel during its operations. On the 
other hand too; supposing that dust breaks off the 
wheel, the dust can not clear the wheel diameter 
because maximum height attained by dust particles 
(from the instantaneous point of break-off) « 1 cm. 


Table 2a: Normal operations at = 1 cm/sec 


Trace 

Gas 

b 

(10“^ gm/sec) 

X 

W„(x) 

W_,(x) 

(drag) 

(sec) 

t^ (vacuum) 
(sec) 

R (drag) 
(cm) 

R (vacuum) 
(cm) 

He 

0.011402 

-0.3262736800 

-0.586591 

-1.573009 

0.008651 

0.010051 

0.003154 

0.005764 

Ar 

0.013177 

-0.3153490738 

-0.542491 

-1.662243 

0.008497 

0.010051 

0.002932 

0.005764 


0.005048 

-0.3578656543 

-0.783090 

-1.253684 

0.009322 

0.010051 

0.004265 

0.005764 


Table 2b: Normal operations at = 1 cm/sec 


Trace Gas 

b 

(10'^ gm/sec) 

tj (drag) 

(sec) 

t^ (vacuum) 

(sec) 

H (drag) 

(cm) 

H (vacuum) 

(cm) 

He 

0.011402 

0.003973 

0.005025473 

0.001504760 

0.002058313 

Ar 

0.013177 

0.003856 

0.005025473 

0.001446236 

0.002058313 


0.005048 

0.004479 

0.005025473 

0.001765332 

0.002058313 


Table 3a: Normal operations at 10 cm/sec. 


Trace 

Gas 

b 

(10“^ gm/sec) 

X 

W„(x) 

W_,(x) 

tj (drag) 
(sec) 

t^ (vacuum) 
(sec) 

R (drag) 
(cm) 

R (vacuum) 
(cm) 

He 

0.011402 

-0.0080385746 

-0.008104 

-6.730092 

0.058954 

0.100509 

0.050243 

0.576498 

Ar 

0.013177 

-0.0037300475 

-0.003744 

-7.622429 

0.057815 

0.100509 

0.043504 

0.576498 


0.005048 

-0.1029401015 

-0.115549 

-3.536842 

0.067776 

0.100509 

0.109913 

0.576498 


Table 3b: Normal operations at 10 cm/sec 


Trace Gas 

b 

(10“^ gm/sec) 

tj (drag) 

(sec) 

t^ (vacuum) 

(sec) 

H (drag) 

(cm) 

H (vacuum) 
(cm) 

He 

0.011402 

0.016721 

0.050254727 

0.047937975 

0.20583131 

Ar 

0.013177 

0.015413 

0.050254727 

0.043096864 

0.20583131 


0.005048 

0.025025 

0.050254727 

0.081468601 

0.20583131 
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Table 4a: Normal operations at 20 cm/sec 


Trace 

Gas 

b 

(10“^ gm/sec) 

x 

W„(x) 

W_,(x) 

tj (drag) 
(sec) 

t^ (vacuum) 
(sec) 

R (drag) 
(cm) 

R (vacuum) 
(cm) 

He 

0.011402 

-0.0000483209 

-0.000048 

-12.46018 

0.109279 

0.201018 

0.100608 

2.305994 

Ar 

0.013177 

-0.0000092725 

-0.000009 

-14.24485 

0.108098 

0.201018 

0.087052 

2.305994 


0.005048 

-0.0139857096 

-0.014186 

-6.073685 

0.120318 

0.201018 

0.227248 

2.305994 


Table 4b: Normal operations at 20 cm/sec 


Trace Gas 

b 

(10'^ gm/sec) 

tj (drag) 

(sec) 

t^ (vacuum) 

(sec) 

H (drag) 

(cm) 

H (vacuum) 

(cm) 

He 

0.011402 

0.035736 

0.100509453 

0.112057611 

0.82332524 

Ar 

0.013177 

0.020158 

0.100509453 

0.099389309 

0.82332524 


0.005048 

0.035736 

0.100509453 

0.209153406 

0.82332524 


8. CONCLUSIONS 

Extremely fine lunar dust clings tenaciously to surfaces 
making the removal efforts quite difficult. Deposition 
of such dust on sensitive surfaces like solar panel, 
optical lenses, etc. causes serious problems. Thus, 
complete understanding of the dynamics of dust during 
rover operations is necessary. However, the non-zero 
gravity of the Moon is manifested as a non-vacuum 
atmosphere near the lunar surface consisting of trace 
gases making determination of trajectory of dust 
difficult. A proposition is arrived to define the 
equations which best describe these dynamics by 
calling for a linear drag force in the equations of 
motion which subsequently leads to transcendental 
equations for which exact closed-form solutions can 
be provided with Lambert fk-function only. However, 
worst case conditions are met considering as if Moon 
does have an ideal vacuum atmosphere. An analysis 
taking Chandrayaan-2 rover is presented and it shows 
that the effect of probable resistance offered by near¬ 
surface lunar atmosphere consisting of an equilibrium 
trace gases is unnoticeable for rovers operating at low 
speed and the effect appears as the rover speed 
increases as expected. It is envisaged that the proposed 
model may also be used for Mars rovers as first order 
approximation. 
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